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Recent observations support the suggestion that short-duration gamma-ray 
bursts are produced by compact star mergers. The X-ray flares discovered 
in two short gamma-ray bursts last much longer than the previously pro- 
posed postmerger energy release time scales. Here we show that they can be 
produced by differentially rotating, millisecond pulsars after the mergers of 
binary neutron stars. The differential rotation leads to windup of interior 
poloidal magnetic fields and the resulting toroidal fields are strong enough to 
float up and break through the stellar surface. Magnetic reconnection-driven 
explosive events then occur, leading to multiple X-ray flares minutes after the 
original gamma-ray burst. 

Gamma-ray bursts (GRBs) are flashes of gamma rays occurring at the cosmological dis- 
tances. They fall into two classes (i): short-duration (< 2 s) hard-spectrum bursts and long- 
duration soft-spectrum bursts. Long GRBs result from core collapses of massive stars (2), and 
short GRBs appear to be produced in mergers of neutron star binaries or black hole-neutron star 
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binaries (5-9). Recently thanks to accurate localizations of several short GRBs (3,6,8) by Swift 
and High Energy Transient Explorer-2 (HETE-2), the multi-wavelength afterglows from these 
events have been detected and the associated host galaxies have been identified. The observa- 
tions provide a few pieces of evidence in favor of the binary compact object merger origin of 
short GRBs (10-12). Because it takes ~ 0.1 — 1 billions years of gravitational wave radiation 
before the binary coalesces, at least some short ORB host galaxies should contain a relatively 
old stellar population. Because neutron stars in the binary system usually receive a very high 
natal velocity, the merger site is preferably at the outskirt of the host galaxy, and the circumburst 
medium density is likely low. These characteristics have been revealed by recent observations: 
First, the identified elliptical galaxies associated with GRB 050509B (3,4) and ORB 050724 
(8,9) suggest that these hosts are early type galaxies with a low star-formation rate, ruling out 
progenitor models invoking active star formation. Second, the nondetection of any supernova 
signal from GRB 050709 indicates that short bursts are not associated with collapses of mas- 
sive stars (5, 7). Third, afterglow modeling of GRB 050709 suggests a low density environment 
(13), which is consistent with that of the outskirt of the host galaxy or that of an intergalactic 
medium. 

However, the above merger origin was recently challenged by the discovery of X-ray flares 
occurring after two short bursts. X-ray flares were discovered to occur at least ~ 100 s after the 
triggers of the short GRB 050709 (J) and GRB 050724 (8). These flares require that the central 
engine is in long-lasting activity. This requirement conflicts with the current models involving 
neutron star-neutron star mergers (14,15) or neutron star-black hole mergers (16), because all 
these models are attached to a common postmerger picture that invokes a black hole surrounded 
by a torus. The predicted typical time scales for energy release are much shorter than > 100 s 
as observed in GRBs 050709 and 050724. Therefore, understanding the origin of X-ray flares 
from short bursts is currently of great interest. Here we show that such flares can be produced 
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by differentially rotating, millisecond pulsars with typical surface magnetic fields that occur 
after the mergers of binary neutron stars. 

In the conventional scenarios of short bursts (10-12), after the merger of a neutron star 
binary, a stellar-mass black hole is formed with a transient torus of mass ~ 1 — 10% of the 
total. These scenarios are valid if the total mass (~ 2.5 — 2.8Mq, where Mq is the solar 
mass) of the postmerger object is larger than the maximum mass of a nonrotating Tolman- 
Oppenheimer-Volkoff neutron star, Mmax.o- This is valid if the nuclear equation of state (EOS) 
is soft to moderately stiff (17). However, the total mass of the postmerger object is smaller than 
Mmax.o for very stiff EOSs (e.g., as predicted by mean field theory) (17). Timing observations 
of the millisecond pulsar J0751+1807 in a circular binary system with a helium white-dwarf 
companion (18) reveal the existence of a neutron star with mass of 2.1 ± 0.2Mq (at the la 
confidence level). This measurement implies that the maximum mass of nonrotating neutron 
stars must be larger than 2.1Mo so that stiff EOSs are favored. Furthermore, recent general 
relativistic numerical simulations (17,19) have shown that for stiff to very stiff nuclear EOSs, 
the postmerger object is indeed a differentially rotating massive neutron star with period of 
~ 1 ms, because uniform rotation and differential rotation can support a maximum mass ~ 
20% and ~ 50% higher than Mjnax,o> respectively. It is therefore reasonable to assume the 
existence of a differentially rotating millisecond pulsar after a double neutron star merger. Such 
a pulsar should also be surrounded by a hot torus with mass ~ 0.01 — O.IMq. Similar to the 
previous scenarios, a short burst may be produced by the Parker instability in the torus (11) or 
the annihilation of neutrinos emitted from the torus (12). 

After the GRB trigger, differential rotation starts to wind the interior magnetic field into a 
toroidal field (20,21). To represent physical processes of windup and floating of the magnetic 
field, we consider a simple two-component model in which the star is divided into two zones 
with a boundary at the radius i?c — 0.5 i?* (where i?* is the stellar radius): the core and the shell 
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components. Their moments of inertia are Ic and Is and their angular (rotation) velocities are 
Qc and Qg, respectively. The differential angular velocity is then AO = Oc — and its initial 
value (marked by a subscript zero) is taken as (A0)o = ^o^s,o (where Aq is the ratio of the 
initial differential angular velocity to the shell's initial angular velocity). If the radial magnetic 

field component is Bj., then the toroidal field component increases as 

^ = (AQ)S,. (1) 

There is a magnetic torque , — {2/3)RlBj.B^, acting between the core and shell (22). This 
torque opposes the differential rotation. Another torque results from magnetic dipole radiation, 
Td = 2B'^Rlill/ (3c^), where c is the speed of light and Bg = eB^ (here e is defined by the ratio 
of the effective surface dipole field strength to the radial field strength). Under action of these 
two torques, the angular velocities of the shell and the core components evolve according to 

Is^^Trr^-Ta (2) 

and 

dVtc 

— Tjfi, (3) 

respectively. The torque from magnetic dipole radiation can be neglected if B^ ^ i?,.e^ (i?* / Rc)^ 
X {R^Q,s/c)^. This condition is easily satisfied at the time t to = e^74Q^0~Q (i?*/i?c)^ (where 
to is ~0.2 ms for typical parameters). Thus, from equations (1-3), we obtained 



dt^ SIJs 

where / = /c + is the total moment of inertia of the star. Letting 



RtB^m): (4) 
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-RlB^ ~ 2.3 X 10^(e/0.3)S-8^ s, (5) 



.34/, 

where 7c — 4 = 10^^ g cm^ and i?* = 10^ cm are taken and § is in units of 10^ G, we found 
a solution of equation (4), 

An = Aofls,QCOs{t/T). (6) 
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This indicates that differential rotation would behave as a resonator if there is no energy dissi- 
pation. 

The increasing toroidal field becomes unstable because of the buoyancy effect when — 
Bjj ~ 10^^ G (20). This corresponds to the time, 

= o f o - X 10\e/0.3)B-lA^'P,,o,u.s, (7) 

where Ps fl,ms is the initial spin period of the shell component in units of milliseconds. Compar- 
ing equations (5) and (7), we see that tb is substantially less than r for ^Ps,o,ms < 1 and that 
the differential angular velocity AQ is approximately constant until the time ti,. At this time, 
the buoyant force is just equal to the force from antibuoyant stratification existing in the star. As 
the time increases, the buoyant force acting on the toroid would begin to exceed the antibuoyant 
force and the toroid will float up toward the stellar surface. The net force density acting on this 
toroid is given by 

BrBbAoQs^o{t — tb) 
f» 9^ (8) 

where Cg is the speed of sound of the embedding medium and g is the surface gravity. In terms 

of equation (8) and Newton's second law, we obtained the buoyancy timescale for the toroid to 

float up and penetrate through the stellar surface: 
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BrBbgAons,o\ y I ) s,8 .,0,ms , V J 

where p ~ 10^^ g cm~^ is the mass density of the embedding medium, and the typical values of 
the speed of sound and the surface gravity are 10^°cms~^ and lO^^cms"^, respectively. This 
timescale is much shorter than tb, suggesting that the toroid, after its field strength reaches Bb, 
would rapidly float up to the stellar surface. 

Once penetrating through the surface, the toroidal fields with different polarity may recon- 
nect (20), giving rise to an explosive event. Its energy is 

Eb = ^Vb ~ 1.6 X 10^1 ergs , (10) 
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where Vb and K are the toroid's volume and stellar volume, respectively. This energy depends 
on the toroid's volume rather than on the initial magnetic field and the stellar spin period. An 
upper limit to the outflow mass ejected is estimated by 

Mb,r.^ = m/g ~ 0.9 X lQ-'MQ{e/Q.3)-'/'B%'ATp;^lt . (11) 

Because of an initial huge optical depth, the outflow will expand relativistically and its minimum 
average Lorentz factor is 

r(.,min ^ 1.0 X 104(e/0.3)^/35-|/'Ao'/'P,'/o!ms- (12) 

The X-ray flares observed at tflare ~ ^6 ^ 100 s after GRBs 050709 and 050724 require that 
the surface magnetic field of a central pulsar i?s ~ 4.8xlO^°(e/0.3)ylo ^Ps^o,ms(^flare/100s)~^ G. 
For typical values (19,22) of the model parameters (i.e., Aq ~ 1, Ps,o ~ 1 ms, and e ~ 0.3), this 
field strength is in the range of the surface magnetic fields of isolated pulsars. Furthermore, it is 
characteristic of the stellar magnetic field that has decayed in ~ 0.1 — 1 billion years before the 
merger of a neutron star binary (25). Inserting this field into equation (12), we found the min- 
imum average Lorentz factor of the outflow from a magnetic-reconnection-driven explosion, 
rb.min ~ 160(tflare/100s)^/^, showing that the outflow is ultrarelativistic. After the end of this 
event, a similar windup of the interior magnetic field with Br would start again following the 
same processes described above, leading to another explosion. 

Collisions among the outflows with different Lorentz factors would produce late internal 
shocks and X-ray flares (24,25). These shocks must produce lower energy photons than did the 
earlier internal shocks during the prompt GRB phase. For the internal shock model, the charac- 
teristic synchrotron frequency is Pm oc L^^'^R'^ oc L^/^F^^^t"^ (where L is the luminosity, i?sh 
is the shock radius, F^ is the bulk Lorentz factor, and 5t is the time interval between two adjacent 
energy shells that the central engine ejects). The late, soft flare is the result of the combination 
of a lower luminosity and a longer time interval (than that of the prompt emission, where 5t ~ 
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in our flare model). In addition, as the stellar differential rotation weakens (i.e., Aq decreases), 
the time interval 6t and the outflow's Lorentz factor increase (see equations 7 and 12). Be- 
cause the maximum flux density of the synchrotron radiation scales as -Fi, max oc F^^ (24), the 
flux density at frequency u is ^ F^,ma^ii^/i^m)~^^~^^^'^ oc F'^'^+p^ 5^^'^^'^ for i/ > i/„ in 
the slow-cooling case (where p is the spectral index of the shock- accelerated electrons) (26). 
Thus, the flare occurring at later times has a smaller flux density because of the larger Lorentz 
factor and longer time interval. This result is consistent with the observed reduced flaring activ- 
ity of GRB 050724. Therefore, our model can provide a self-consistent explanation for all the 
observations including the energetics (see equation 10) and the temporal and spectral properties 
of the X-ray flares. 

Generally speaking, the surface magnetic field of the postmerger pulsar could have a wider 
range than the preferred value invoked here to interpret the ~ 100 s flares in GRBs 050709 and 
050724. For stronger fields, this would give rise to multi-peaks in the prompt phase (as observed 
in some short GRBs) or, if the flares are not bright enough, they may be masked by the steep 
decay component of the prompt emission tail (25). For weaker fields, the putative flares occur 
much later and are energetically insignificant. This would give rise to smoother X-ray afterglow 
lightcurves as observed in several GRBs observed by Swift (e.g. GRB 050509B) (3). 

X-ray flares were observed in nearly a half of long Swift bursts (27,28). Even though the 
two classes of bursts have different progenitors (namely coUapsars for long bursts and binary 
neutron star mergers for short bursts), similar temporal properties (e.g. peak times and temporal 
indices before and after the peaks) suggest that the X-ray flares may have a common origin. 
Therefore, we suggest that some long bursts may originate from moderately-magnetized mil- 
lisecond pulsars with hyperaccreting accretion disks after the collapses of massive stars and 
their X-ray flares are the result of strong interior differential rotation of these pulsars. The 
differences in duration, energetics and spectrum for the two classes of bursts would be due to 
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different accretion disks, e.g., a transient torus for short bursts (10-12) and a fall-back accretion 
disk for long bursts (29,30). When the surface magnetic fields are strong enough, the spin down 
of this central engine pulsar would provide energy injection to the postburst relativistic outflow 
(31), which could interpret the late X-ray humps detected in many GRBs (25,28). 
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